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ABSTRACT 

A robust connection between the drag on surface-layer winds and the stratospheric circulation is demon- 
strated in NASA’s Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM). Specifically, 
an updated parameterization of roughness at the air-sea interface, in which surface roughness is increased for 
moderate wind speeds (4-20 ms -1 ), leads to a decrease in model biases in Southern Hemispheric ozone, 
polar cap temperature, stationary wave heat flux, and springtime vortex breakup. A dynamical mechanism is 
proposed whereby increased surface roughness leads to improved stationary waves. Increased surface 
roughness leads to anomalous eddy momentum flux convergence primarily in the Indian Ocean sector (where 
eddies are strongest climatologically) in September and October. The localization of the eddy momentum flux 
convergence anomaly in the Indian Ocean sector leads to a zonally asymmetric reduction in zonal wind and, 
by geostrophy, to a wavenumber-1 stationary wave pattern. This tropospheric stationary wave pattern leads to 
enhanced upward wave activity entering the stratosphere. The net effect is an improved Southern Hemi- 
sphere vortex: the vortex breaks up earlier in spring (i.e., the spring late-breakup bias is partially ameliorated) 
yet is no weaker in midwinter. More than half of the stratospheric biases appear to be related to the surface 
wind speed biases. As many other chemistry-climate models use a similar scheme for their surface-layer 
momentum exchange and have similar biases in the stratosphere, the authors expect that results from 
GEOSCCM may be relevant for other climate models. 


1. Introduction 

The interaction between the surface and the lowest 
levels of the atmosphere is a crucial process in nature 
and in any climate model. The midlatitude surface 
westerlies are maintained against surface drag by eddy 
momentum fluxes (Held 1975; Edmon et al. 1980; Vallis 
2006, section 12.1). Changes in surface drag can cause 
meridional shifts of the midlatitude, eddy-driven jet 
(Robinson 1997; Chen et al. 2007; Kidston and Vallis 
2012). The flux of enthalpy into the atmosphere at the 


Corresponding author address: Chaim I. Garfinkel, Department 
of Earth and Planetary Science, The Johns Hopkins University, 
Baltimore, MD 21218. 

E-mail: cig4@jhu.edu 

DOI: 10.1175/JAS-D-12-0242.1 


bottom of a hurricane regulates subsequent growth of 
the hurricane (Molod et al. 2013). Increased surface drag 
leads to more efficient eddy momentum fluxes into 
the jet core, so that even as the increased roughness is 
reducing winds, eddies try to counter this effect by more 
efficiently transporting momentum into the jet (Barnes 
and Garfinkel 2012). 

The exchange of momentum, moisture, and sensible 
heat between the ocean surface and the atmosphere 
occurs on spatial and temporal scales far finer than any 
atmospheric general circulation model (GCM) can di- 
rectly simulate. Many models therefore rely on Monin- 
Obukhov similarity theory (MOST) to specify air-sea 
momentum exchange as a function of bulk winds, tem- 
perature, and humidity. However, the coefficients un- 
derlying MOST, and specifically those relating the 
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air-sea roughness to the turbulence in the near-surface 
flow, are observationally derived. Garfinkel et al. (2011) 
found that updating these coefficients based on the most 
recently available observational data for wind speeds 
between 4 and 20 m s _1 leads to improved surface winds 
over the Southern Ocean in the National Aeronautics 
and Space Administration (NASA)’s Goddard Earth Ob- 
serving System Chemistry-Climate Model (GEOSCCM). 
In addition, they found that the updated parameteriza- 
tion leads to improved zonal surface stress on the ocean 
and eddy momentum flux convergence aloft. They also 
mentioned, but did not discuss in detail, the improve- 
ment of the Southern Hemispheric (SH) stratospheric 
circulation. Here, we expand on Garfinkel et al. (2011) 
and discuss these improvements in the SH stratosphere. 

A persistent bias in most chemistry-climate models 
(CCMs) is the late breakup of the SH stratospheric polar 
vortex in spring [September-November; cf. chapter 4 of 
SP ARC-CCM V al (2010); Butchart et al. (2011)]. This 
bias tempers the conclusions that can be drawn from 
model output data. For example, details regarding the 
timing of the tropospheric impacts of ozone depletion 
(e.g., Arblaster and Meehl 2006; Polvani et al. 2011) and 
the modeled ozone trends in the SH lower stratosphere 
(e.g., Stolarski et al. 2006) are likely also biased in many 
models. In the GEOSCCM, the springtime vortex bias 
has been present in multiple model generations (e.g., 
Stolarski et al. 2006; Hurwitz et al. 2010) despite im- 
provements to model parameterizations, including the 
gravity wave parameterization. Hurwitz et al. (2010) 
showed that the delayed breakup in GEOSCCM is re- 
lated to overly weak SH lower-stratospheric heat flux 
in spring, and in particular, to weak stationary planetary 
waves in the troposphere. 

In this paper, we show that a by-product of the im- 
proved air-sea roughness parameterization is improved 
tropospheric stationary waves. Once the tropospheric 
stationary waves are improved, stratospheric sensible 
heat flux, and in particular that associated with station- 
ary planetary waves, is improved as well. The polar 
vortex breaks down earlier in spring, and biases in 
stratospheric ozone are partially ameliorated. As many 
other atmospheric CCMs appear to implement a similar 
air-sea roughness parameterization for momentum ex- 
change as GEOSCCM and base it on similarly old and 
scant empirical data (e.g., Garfinkel et al. 2011, their 
Table 1), we expect that the reduction in model bias 
shown here may be common to other CCMs as well. 

This paper is organized as follows. After explaining 
the methodology in section 2, we document the improve- 
ments in stationary waves, heat flux, polar cap temper- 
ature, 60°S winds, and stratospheric ozone in sections 3 
and 4. We then present a mechanism whereby the change 


in air-sea roughness over an essentially zonally sym- 
metric ocean can lead to improved tropospheric station- 
ary waves in section 5. Finally, section 6 demonstrates 
that more than half of the stratospheric biases appear to 
be related to the surface wind biases. 

2. Methods and data 

The GEOSCCM experiments performed use the 
Goddard Earth Observing System, version 5 (GEOS-5), 
atmospheric GCM (Rienecker et al. 2008) coupled to a 
comprehensive stratospheric chemistry module (Pawson 
et al. 2008). The model has 72 vertical layers, with a 
model top at 0.01 hPa, and all simulations discussed here 
were performed at 2° latitude X 2.5° longitude hori- 
zontal resolution. A 25-yr run at 1° X 1.25° horizontal 
resolution without interactive stratospheric chemistry 
has been performed to test sensitivity to resolution, and 
results are qualitatively similar (not shown). SPARC- 
CCMVal (2010) grades highly the representation of 
the Southern Hemisphere stratosphere by an earlier 
version of GEOSCCM as compared to observations and 
to the multimodel mean of an ensemble of CCMs. 

Details of the surface-layer scheme used in GEOSCCM 
can be found in the appendix of Helfand and Schubert 
(1995), and a description of the update of the air-sea 
roughness parameterization can be found in Garfinkel 
et al. (2011). Briefly, Garfinkel et al. (2011) showed 
that the local effect of updating the parameterization 
was 1) to increase the roughness and the momentum ex- 
change coefficient and 2) to decrease surface wind speed, 
over the oceanic regions where the wind speed exceeds 
4ms" 1 . Here, the effect of updating the air-sea roughness 
parameterization on the stratospheric circulation is eval- 
uated in the GEOSCCM. The air-sea roughness pa- 
rameterization in the CONTROL (i.e., lower roughness) 
simulation is different from that in the UPDATED 
[i.e., updated as in Garfinkel et al. (2011) with higher 
roughness] simulation; all other boundary conditions 
and model settings are identical in the two simulations. 
Three 30-yr “CONTROL” and 30-yr “UPDATED” 
model integrations were performed. While the three 
ensemble members differ in their gravity wave drag 
scheme [one follows Rienecker et al. (2008), while two 
follow Molod et al. (2012)], the key results presented in 
this paper are insensitive to the choice of gravity wave 
drag scheme. Thus, in the discussion of our results, we do 
not distinguish between the three ensemble members. 

Greenhouse gas and ozone-depleting substance con- 
centrations represent the year 2005 in all experiments 
discussed here. The sea surface temperatures in the 
years 1980-2010 force each 30-yr pair. Variability re- 
lated to the solar cycle and volcanic aerosols are not 
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considered. The parameterization of enthalpy fluxes at 
the air-sea interface is held fixed between UPDATED 
and CONTROL. The key point is that these simulations 
isolate the impact of updating the air-sea roughness 
parameterization for momentum. 

a. Validation 

NASA’s Modern-Era Retrospective Analysis for Re- 
search and Applications (MERRA; Rienecker et al. 
2011) is used to validate the model’s height, eddy fluxes, 
and zonal wind away from the surface; results are nearly 
identical when we use 40-yr European Centre for 
Medium-Range Weather Forecasts (ECMWF) Re- 
Analysis (ERA-40; Uppala et al. 2005). The model’s 
ozone is compared to the patched and interpolated 
satellite-based product available online (at http://www. 
bodekerscientific.com/data/total-column-ozone) [updated 
from Bodeker et al. (2001, 2005)]. Modeled surface 
winds and surface stress are validated against version 2 
of the Goddard Satellite-Based Surface Turbulent Fluxes 
(GSSTF) data [available at http://disc.sci.gsfc.nasa.gov/ 
precipitation/data-holdings/access/gsstf2.0.shtml and 
described in Chou et al. (2003)]. For ozone and quan- 
tities derived from MERRA, we begin the validation 
period in 1989 since the model is forced with 2005 
ozone-depleting substances and we wish to restrict our 
analysis to observed years with comparable levels of 
ozone-depleting substances. 

b. Diagnostics 

In section 5, we will use a vorticity budget to diagnose 
the role of eddies in establishing a tropospheric zonal 
wavenumber-1 (hereafter wave-1) response to the up- 
dated air-sea drag. Here, we introduce the budget used 
to demonstrate this effect. Barnes and Hartmann (2010a, 
hereafter BH10) and Barnes and Hartmann (2010b) 
show that the interactions between temporal eddies 
and the time mean flow in zonally confined regions can 
be diagnosed quantitatively by using the vorticity budget 
(Holton 2004): 

J t = [-(£ + /)V • U - A «] stretching + [-v • (u£)] eddy 

+ [-u-V(f+/)-u-Ve] wave 

+ {-V.[u£ + /)]} clim +^. 

( 1 ) 

The overbar represents the time mean and the caret 
represents the deviation from the time mean, f is the 
vorticity, /is the Coriolis parameter, u is the vector wind, 
and T is the drag as applied to vorticity and is pro- 
portional to the surface stress [e.g., section 7 of Barnes 
and Garfinkel (2012)]. BH10 demonstrate that when the 


vorticity forcing terms in upper levels project onto vor- 
ticity anomalies in lower levels, lower-level vorticity 
anomalies can be maintained against damping. This 
budget was also used to demonstrate that eddies re- 
spond to an external forcing in Garfinkel and Hartmann 
(2011). See BH10 for more details. 

The vorticity budget is used here to diagnose the role 
of tropospheric eddies in the tropospheric response to 
the updated surface drag. The terms in the vorticity bud- 
get are computed for the CONTROL and UPDATED 
runs. The difference in the forcing terms between the 
UPDATED and the CONTROL runs represents the 
generation of tropospheric vorticity anomalies in re- 
sponse to increased surface drag. Specifically, we will 
show that a wave-1 surface stress anomaly T is balanced 
by the other terms in the budget. 

A two-tailed Student’s t test is used to assess statistical 
significance where indicated. Each year is taken as 1 
degree of freedom, except when we discuss the high- 
frequency poleward eddy momentum flux where each 
month is taken as 1 degree of fre edom. (Quasi-) sta- 
tionary heat flux is computed from Vso^so^ where v$ 0 is 
the deviation of the meridional wind from the zonal 
mean averaged over September and October (SO), T$ 0 
is like u| Q , but for temperature, and the overbar denotes 
a zonal average. Heat flux is used as a proxy for wave 
activity flux. High-frequency eddy momentum fluxes are 
computed by applying a 7-day high-pass ninth-order But- 
terworth filter applied separately to the meridional and 
zonal velocities before computing their product. 

We now address the impact of the change in the air- 
sea roughness parameterization on bulk quantities in the 
model. 

3. Improvements in the zonally averaged 
circulation 

We first focus on the bias in polar cap temperature 
(75°S and poleward) of the CONTROL climate and 
its improvement upon updating the air-sea roughness 
parameterization. The climatological polar cap temper- 
ature in observations and in the model are shown in 
Figs, la and lb, and the difference between the two cli- 
matologies is shown in Fig. lc. In the CONTROL cli- 
mate, modeled temperatures are too warm in autumn and 
winter yet are too cold in spring. When we update the air- 
sea roughness parameterization, we find that while the 
bias during winter is unchanged, the biases during fall and 
spring are reduced (Fig. Id). In the lower stratosphere, the 
springtime bias is reduced by around one-third. 

Subpolar (60°S) zonal-mean zonal wind is also improved 
when we update the air-sea roughness. The climatological 
subpolar zonal winds in observations and in the model 
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Fig. 1. Polar cap temperature (area- weighted average, 75°S and poleward) evolution over the course of the year in 
(a) CONTROL, (b) reanalysis data, (c) CONTROL minus reanalysis, and (d) UPDATED minus CONTROL. The contour 
interval is 5 K in (a),(b), 2 K in (c), and 0.5 K in (d). Regions with anomalies significant at the 95% (99%) level are orange for 
positive anomalies or light blue for negative anomalies (red or dark blue, respectively) in (d). The zero contour is omitted. 


are shown in Figs. 2a and 2b and the difference between 
the two climatologies is shown in Fig. 2c. In the CON- 
TROL climate, winds are too strong in spring and 
summer. When we update the air-sea roughness pa- 
rameterization, we find that westerlies weaken in all 
seasons in the troposphere and lower stratosphere, but 
the effect is strongest in SH spring. In the middle and 
upper stratosphere, winds are decreased only in spring 
and summer. The biggest easterly anomalies occur where 
the CONTROL experiment is biased toward too-strong 
westerlies (in late spring and early summer). Note that 
the midwinter vortex is not significantly affected by the 
change in the air-sea roughness parameterization, but the 
vortex breaks up earlier in spring. The earlier vortex 
breakup is clearly shown in the timing of the springtime 
transition to easterlies in the SH at 60°S, one of the 
metrics identified by SPARC-CCMVal (2010). Figure 3b 
shows that by this metric as well, GEOSCCM is now 
improved by approximately one week near lOhPa (33% 
of the bias), though model biases do still exist. 1 


1 The magnitude of this bias is much smaller at 1° X 1.25° hori- 
zontal resolution (not shown). 


Zonally averaged polar cap (75°S and poleward) col- 
umn ozone is also improved in GEOSCCM in late spring 
(Fig. 3c). In November and December, there is signifi- 
cantly more ozone in the SH polar region. More than 
half of the bias in GEOSCCM in these months is elim- 
inated by the updated air-sea roughness parameteriza- 
tion, though the large variability in the observational 
data complicates our interpretation. The improvement 
in spring column ozone is consistent with the earlier 
vortex breakup. The high ozone bias of the model (e.g., 
Pawson et al. 2008) complicates the interpretation of the 
results in January. The change in May and June appears 
to be related to a change in transport across the polar 
vortex: column ozone is significantly increased between 
30° and 60°S (not shown). There is no net effect if we 
average from 30°S to the pole. A detailed investigation 
of the change in late fall is beyond the scope of this work. 

Interannual variability of the spring vortex breakup is 
also improved. To demonstrate this, we show the stan- 
dard deviation of polar cap column ozone among the 90 
years in CONTROL and in UPDATED in Fig. 3d. 
There is too little interannual variability in CONTROL. 
In other words, the breakup date of the ozone hole is too 
tightly clustered around the same time period (which in 
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Fig. 2. As in Fig. 1, but for 60°S zonal wind evolution over the course of the year. The contour interval is 8 m s 1 in 

(a),(b), 2 ms -1 in (c), and 1ms -1 in (d). 


turn is also biased as shown in Fig. 3c). In contrast, more 
interannual variability is present in spring in UPDATED, 
though the springtime peak in variability is still delayed 
as compared to observations. In summary, the zonally 
averaged circulation, and in particular polar cap temper- 
ature, 60°S zonal wind, and column ozone, is improved in 
midspring through early summer by updating the air-sea 
roughness parameterization. 

4. Improvements in the zonally asymmetric 
circulation 

We now discuss zonally asymmetric improvements in 
September-October. We focus on September-October 
because it precedes the largest improvements in the 
zonally averaged stratospheric circulation (Hurwitz 
et al. 2010). Figures 4 and 5 show that surface wind speed, 
upper-tropospheric zonal wind, upper-tropospheric height, 
and column ozone are all improved by the updated pa- 
rameterization in September-October. As discussed in 
Garfinkel et al. (2011), wind speed at the surface is im- 
proved. This can be seen by noting that the regions of 
largest bias in the CONTROL run in Fig. 4b have op- 
posite changes in Fig. 4c. In addition, the pattern of the 
change resembles zonal wave 1: winds are decelerated 


most strongly in the Pacific sector (Fig. 4c). Aloft, the 
bias in the model, and the improvement upon updating 
the air-sea roughness parameterization, also resemble 
a strong wave-1 pattern (Figs. 4d-f). Associated with the 
reduction in zonal wind over the Southern Ocean is 
a subpolar high and subtropical low (Figs. 5a,b) by ge- 
ostrophy. Since the reduction in zonal wind is concen- 
trated in the Pacific sector, the anomalous ridge is 
strongest in this region. 

To show the vertical structure of the height anoma- 
lies, we compute the deviation of height from the zonal 
mean (i.e., the zonally asymmetric component), weight 
it by the square root of the air density [as in Matsuno 
(1970) and in section 12.3.1 of Holton (2004)], and 
then plot the result in Fig. 6. Stationary planetary (in 
particular wave-1) waves are too weak in CONTROL 
(Figs. 6b,c,e), but the bias is partially ameliorated in 
UPDATED throughout the troposphere and strato- 
sphere (cf. the color hues in Figs. 6d and 6e). The im- 
provement generally resembles wave-1; consistent with 
Figs. 5a and 5b, positive height anomalies are larger in the 
Pacific sector in the troposphere (Fig. 6f). In the strato- 
sphere, the anomaly has westward tilt with height, which 
is indicative of enhanced wave activity flux into the 
stratosphere. 
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Fig. 3. Evolution of (a) polar cap temperature at lOOhPa, (b) the date on which 60°S zonal wind reaches 0ms -1 
(i.e., the vortex breakup), (c) polar cap column ozone, and (d) variability in polar cap column ozone, in the 
UPDATED, CONTROL, and reanalysis data. The range of observational variability is indicated with gray shading 
or error bars. Days in which UPDATED and CONTROL are significantly different from each other at the 95% level 
are designated with a star below the curves or to the left of the curves. The polar cap average is defined as the area- 
weighted average 75°S and poleward. Monthly averaged data is used for the observations for (c),(d), and daily data 
are used for the observations for (a),(b). For (d), an F test is used to assess significance. 


Zonally asymmetric improvements in ozone are also 
present in September-October. Namely, the bias in col- 
umn ozone, as well as its improvement upon increasing the 
air-sea roughness, also bears a wave-1 pattern (Figs. 5c, d). 
A comparison of the phasing of the ozone improvement 
in Figs. 5c and 5d and the height improvement in Fig. 6f 
suggests that the geostrophic flow associated with the 
height improvement can, by advecting ozone, lead to 
improved ozone as well. Namely, poleward motion 
south of the Indian Ocean advects ozone-rich subpolar 
air poleward, while equatorward motion south of South 
America advects ozone-poor polar air equatorward. 
Even though the zonally averaged improvements in 
ozone do not begin until November (Fig. 2c), zonally 
asymmetric improvements are present earlier. 

Next, we consider the change in heat flux (Fig. 7). In 
the stratosphere, the model has too little heat flux 
(Figs. 7a, b). This bias in the stratosphere is partially 
ameliorated when we update the air-sea roughness 


parameterization (Fig. 7c). 2 Nearly all of the bias in the 
model, as well as the improvement in UPDATED, is 
associated with the stationary waves as opposed to the 
transient waves (Figs. 7d-f). It is not surprising that the 
improvement in stationary wave heat flux and in height 
should be linked. Garfinkel et al. (2010) and Smith and 
Kushner (2012) argue that if an anomalous ridge is col- 
located with a ridge of the climatological planetary wave 
pattern (as is evident in the troposphere in Fig. 6f), then 
tropospheric upward wave activity flux will be enhanced. 
Figure 8 confirms that the change in heat flux in the 
stratosphere is associated almost entirely with wave 1. 
However, the too-low heat flux bias in GEOSCCM dis- 
cussed in Hurwitz et al. (2010) is not removed entirely. In 
the troposphere as well (where synoptic transient heat 


2 The magnitude of this bias is much smaller at 1° X 1.25° hori- 
zontal resolution (not shown). 
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Fig. 4. (a),(d) Control run climatology, (b),(e) model bias, and (c),(f) improvement upon updating the air-sea roughness parameter- 
ization in (a)-(c) surface wind speed, and (d)-(f) 250-hPa zonal wind, during September and October. The contour intervals are 2.25 ms -1 
in (a), 0.7 m s -1 in (b),(c), 9 m s -1 in (d), and 1 m s -1 in (e),(f). Regions with anomalies significant at the 95% (99%) are light (dark) red or 
blue in (e) and (f), with red denoting positive anomalies and blue denoting negative anomalies. The zero contour is omitted and negative 
contours are dashed and thick. 


flux dominates stationary and planetary wave heat flux 
climatologically), planetary wave heat flux is improved 
(Fig. 9). In particular, wave-1 heat flux between 50° and 
70°S is significantly increased (Fig. 9d). We defer a discus- 
sion of the changes for synoptic wavenumbers in the tro- 
posphere and for forecasting skill, and a comparison of our 
results to those of Branscome et al. (1989), for future work. 

In summary, the updated air-sea roughness parame- 
terization leads to a strong wave-1 response in September- 
October. Zonally asymmetric height, zonal wind, and eddy 
heat flux are all improved. The improved heat flux in 


September-October then drives the zonally symmetric 
improvements in the late spring and early summer strato- 
spheric circulation discussed in section 3. 

5. How does surface roughness impact 
stationary waves? 

We now seek to mechanistically explain the origin of 
these zonally asymmetric improvements in the upper- 
tropospheric and stratospheric circulation. In particular, 
we seek to answer the following questions. 



Fig. 5. (a),(c) Model bias and (b),(d) improvement upon updating the air-sea roughness parameterization in (a),(b) 250-hPa height, and 
(c),(d) column ozone, during September and October. The contour intervals are as follows: 20 m in (a),(b), 8 Dobson units (DU) in (c), and 
4 DU in (d). Regions with anomalies significant at the 95% (99%) are light (dark) red or blue in (c),(d), with red denoting positive 
anomalies and blue denoting negative anomalies. The zero contour is omitted and negative contours are dashed and thick. 
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Fig. 6. Longitude vs altitude cross section of the zonally asymmetric component of geopotential height at 60°S during September and 
October in (a) UPDATED, (b) CONTROL, (c) reanalysis data, (d) UPDATED minus reanalysis, (e) CONTROL minus reanalysis, and 
(f) UPDATED minus CONTROL. The contour interval is 15 m for (a)-(c) and 7.5 m for (d)-(f). Height has been weighted by \/p/p s , the 
ratio of the pressure at each level to the surface pressure, as described in the text. Additional panels are shown in order to ease comparison 
with Fig. 7 of Hurwitz et al. (2010). The zero contour is thick. 


(i) What is the mechanism whereby a change in air-sea 
roughness over the essentially zonally symmetric 
Southern Ocean can lead to a wave-1 tropospheric 
response? 

(ii) What controls the phase of the wave-1 pattern? 

One might hypothesize that the response to drag is pro- 
portional to the local climatological wind speed and thus 
bears a slight zonal wave-1 structure as the underlying 
Southern Ocean winds bear a slight wave-1 structure 
(cf. Fig. 4a; such a prediction follows from running the 
surface layer scheme offline from the rest of the model.). 
However, the strongest climatological surface winds are 
in the Indian Ocean sector, while the strongest response is 
in the Pacific sector and not in the Indian Ocean sector 
(Fig. 4a vs Fig. 4c). Thus, a simple scaling argument based 
on the climatological wind speed fails to predict the re- 
sponse. The rest of this section presents a more plausible 
mechanism as to how the change in roughness can create 
a wave-1 anomaly in the troposphere. 

As mentioned in the introduction, Barnes and Garfinkel 
(2012) note that as drag on winds is increased in a baro- 
tropic model, eddy momentum fluxes into the jet core are 


enhanced if the drag is quadratic (or cubic, once one ac- 
counts for the increase in the momentum exchange co- 
efficient with wind speed) in wind speed as in GEOSCCM. 
Even as the increased roughness is reducing winds, eddies 
try to counter this effect by more efficiently transporting 
momentum into the jet. This behavior is also evident in 
Fig. 3 of Kidston and Vallis (2012). 

In the SH, the region of eddy momentum flux is 
zonally asymmetric in September-October, as seen in 
Fig. 10a. In particular, Fig. 10a shows that more high- 
frequency eddy momentum flux convergence (EMFC) 
occurs in the Indian Ocean sector than in the Pacific 
sector [consistent with Barnes and Hartmann (2010b)]. 
The wave-1 zonal asymmetry is especially strong during 
SH winter, but it is present in September-October and 
in the annual average as well. It is also present in other 
diagnostics of storm tracks (e.g., Hoskins and Hodges 
2005, their Figs. 2 and 11; Chang et al. 2013, their Fig. 1). 

The bias in EMFC in CONTROL, and the improve- 
ment upon increasing the surface roughness, also bears 
a zonally asymmetric wave-1 pattern (Figs. 10c, e). In 
the Indian Ocean sector, eddies react strongly to the 
weakening of the jet by enhanced roughness (Fig. lOe). 
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Heat Flux 



latitude 


Fig. 7. Zonally averaged sensible heat flux during September and October in (a),(d) CONTROL, (b),(e) reanalysis data, and (c),(f) 
UPDATED minus CONTROL, (a)-(c) The total heat flux and (d)-(f) stationary waves only are shown. The difference upon increasing 
the roughness is shown in (c), (f). The contour interval is 8Kms _1 in (a),(b), 4Kms -1 in (d),(e), and 0.5 K ms -1 in (c),(f). The zero 
contour is thick in (c),(f). 


In the Pacific sector, however, the torque due to the 
eddies is weak and so the decrease in zonal wind because 
of enhanced roughness is larger there. In other words, 
a simple scaling argument applied to the climatological 
eddies (rather than the climatological surface wind speed) 
can explain the response of eddies to increased surface 
drag: the response of eddies to drag is proportional to the 
local intensity of eddy activity and thus features a pro- 
nounced maxima in the Indian Ocean sector. Because of 
this, the strongest change in zonal wind due to eddies is 
in the Indian Ocean sector. Finally, because eddies act to 
accelerate the jet, they counteract the effect of increased 
friction most strongly in the Indian Ocean sector. 

A close examination of Fig. lOe suggests that the 
EMFC is positive only on the equatorward flank of the 
Indian Ocean sector jet and not on the poleward flank of 
the jet — the change resembles a meridional dipole — and 
so it is not clear whether the net effect in the Indian 
Ocean sector is to accelerate or decelerate the jet. There 
are two points to bear in mind when considering this 
issue. First, the jet itself shifts equatorward in response 
to surface friction [see Fig. 4f, Robinson (1997), and 
Kidston and Vallis (2012)]. The equatorward shift is 
amplified here because the polar stratospheric vortex is 
weakened and a weakened vortex tends to shift tropo- 
spheric jets equatorward [see Fig. 1, and also Garfinkel 
et al. (2012), and the references therein]. We therefore 


should expect an equatorward shift in EMFC. Second, 
we have analyzed the total zonally asymmetric forcing of 
vorticity which was introduced in section 2b, and we find 
that the subsequent interaction of the eddies with the 
mean state 1) is more clearly associated with an accel- 
eration in the Indian Ocean sector poleward of the jet 
core and 2) is stronger than the deceleration because 
of the high-frequency eddies in this region. This effect is 
easily seen in Fig. 11, which shows the upper-tropospheric 
height, vorticity, and total forcing (azonal plus zonal) of 
vorticity in the SH averaged from 50° to 70°S. The total 
forcing of vorticity is defined as the right-hand side of 
Eq. (1) except for J 7 , the direct effect of the frictional 
forcing. The direct effect of increasing roughness is to 
increase (i.e., a positive anomaly of) subpolar vorticity 
and geopotential height because of geostrophy (red and 
black curves in Fig. 11). The eddy forcing counteracts this 
change and reduces vorticity in the zonal average (blue 
curve in Fig. 11), but the primary reduction is in the In- 
dian Ocean sector. The eddy fluxes are weakest clima- 
tologically in the Pacific sector as well (e.g., Fig. 10a). The 
net effect is that the height and vorticity response to in- 
creased roughness is strongest in the Pacific sector where 
the eddy fluxes are weakest. 

The changes in eddy momentum flux convergence 
aloft imply anomalous surface stress, as du*v*/dy must 
balance the removal of surface momentum by friction 
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Fig. 8 As in Fig. 1, but for u* 7* at lOOhPa during September-October as a function of zonal wavenumber. The 
contour interval is 2Kms _1 in (a),(b) and 2/3 K ms -1 in (c),(d). The zero contour is thick. 


for a steady-state surface jet (Held 1975; Vallis 2006, 
section 12.1). Figures 10d,f show that zonal surface stress 
is also improved with the updated surface roughness as 
discussed in Garfinkel et al. (2011). In addition, it is clear 
that the change in surface stress bears a wave-1 signa- 
ture. Furthermore, the phase of the wave-1 pattern is 
different than that for, for example, height or zonal 
wind, but is similar to that of EMFC; namely, the largest 
change in surface stress is in the Atlantic/Indian Ocean 
sector and not in the Pacific Ocean sector. It is re- 
markable that the changes in surface stress and eddy 
momentum flux convergence balance each other nearly 
perfectly even without performing a zonal average (cf. 
Figs. lOe and lOf). Hence, we can explain the differences 
in phasing between the surface stress and eddy flux 
changes (which are largest in the Indian Ocean sector) 
and the zonal wind changes (which are largest in the 
Pacific sector): the surface stress changes most strongly 
where the EMFC change is largest, while the surface 
wind changes most strongly where the EMFC change is 
smallest. 

The schematic in Fig. 12 summarizes the observed 
changes in the troposphere. Climatologically, the drag 
and eddies balance each other such that there are 


westerlies over the Southern Ocean. When we update 
the air-sea roughness parameterization, we increase the 
drag in a nearly zonally symmetric manner (and if any- 
thing, the change is strongest in the Indian Ocean sector 
if the air-sea roughness scheme is run offline). This is 
represented by the orange easterly arrows in both the 
Indian Ocean and Pacific sectors. However, eddies act to 
mitigate this easterly anomaly; this effect is represented 
by dark blue arrows. Eddies are stronger climatologi- 
cally in the Indian Ocean sector, and thus, the dark blue 
arrow is longer in the Indian Ocean sector than in the 
Pacific sector. The net effect is that the deceleration of 
the wind, which is represented with a purple-gray arrow, 
is larger in the Pacific sector than in the Indian Ocean 
sector. 

If we revisit the changes in zonal wind and height in 
Figs. 4a-f and 5a-b, it is clear that the zonal wind re- 
duction at the surface and aloft because of increased 
roughness is strongest in the Pacific sector. Eddies ac- 
celerate the flow throughout the midlatitude Indian 
Ocean sector [as suggested by the simpler experiments 
of Barnes and Garfinkel (2012)], but not in the Pacific 
Ocean sector; hence, increased surface roughness leads 
to a wave-1 response. 
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Sensible heat flux (500hPa) decomposed by wavenumber 
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Fig. 9. As in Fig. 8, but for v*T* at 500 hPa during September-October. The contour interval is 0.4 Km s 1 in 
(a),(b) and 0.1333 Km s -1 in (c),(d). The zero contour is thick. 


Now that we have mechanistically explained the ori- 
gin of the wave-1 pattern for the zonal wind change 
(because the eddies do not compensate as strongly in 
the Pacific sector), the rest of the response follows 
simply. The anomalous ridging is stronger in the Pacific 
sector because of geostrophy. This anomalous ridge 
then affects the stationary waves (Figs. 6-8) and leads 
to increased wave driving of the vortex. The increased 
heat flux entering the stratosphere in September-October 
causes an earlier breakup of the stratospheric vortex. 
The key point is that because eddies are less effective at 
counteracting the effect of enhanced ocean roughness 
in the Pacific sector as compared to the Indian Ocean 
sector, the net response will bear a wave-1 structure. 

This mechanism is consistent with the timing of 
the improvement in stationary waves and in polar cap 
temperature. The midwinter SH vortex appears to be 
too strong to be greatly affected by enhanced tropo- 
spheric resolved waves (Charney and Drazin 1961; Scott 
and Haynes 2002; Taguchi and Yoden 2002). Therefore, 
it is to be expected that the late-breakup bias can be al- 
leviated while not exacerbating the midwinter bias once 
resolved stationary waves are improved. In addition, this 
mechanism would suggest that the fall spinup of the 


vortex would not be affected by the updated air-sea 
roughness parameterization. The mechanism relies on 
the existence of stationary waves in the troposphere, 
and in February- April (when the vortex spins up), the 
tropospheric circulation is more zonally symmetric 
[not shown, but similar to Fig. 2 of Hoskins and Hodges 
(2005)]. 

Note that this mechanism appears disconnected from 
the barotropic governor of James and Gray (1986) and 
James (1987). Wave-1 stationary waves are the cause of 
the increased wave driving of the vortex and not higher- 
frequency synoptic waves as might be expected from an 
argument involving baroclinic instability. 

In summary, Barnes and Garfinkel (2012) show that 
eddies can change so as to counteract the effects of in- 
creased friction on the mean flow; that is, eddies act to 
accelerate the jet. The mechanism we present takes this 
just one step further: eddies are effective at compen- 
sating for the increased surface roughness mainly where 
high-frequency eddies are strongest climatologically 
(i.e., in the storm tracks). Because the storm tracks have 
a wave-1 pattern in September-October in the SH, the 
response to increased surface roughness will have a 
wave-1 structure as well. 
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run, (c) CONTROL minus observations, and (e) UPDATED minus CONTROL. Contour interval is 0.75 ms -1 day -1 in (a), 0.3 ms -1 day -1 
in (c), and 0.15 m s -1 day -1 in (e). (b),(d),(f) As in (a),(c),(e), but for eastward surface stress. Contour interval is 5 X 10 -2 N m -2 in (b), 
2 X 10 -2 Nm -2 in (d), and 10 -2 Nm -2 in (f). Regions with anomalies whose statistical significance exceeds 95% are in color in (e),(f). 
The zero contour is omitted and negative contours are dashed. For this figure only, we show the annual average, as the September- 
October average is very noisy (though qualitatively similar) for these quantities. 


6. Stratospheric improvement upon removing 
entirely the surface wind bias 

In all previous sections of this paper, we have con- 
sidered the changes in the model upon updating the air- 
sea roughness coefficients to match recent observations 
(e.g., Garfinkel et al. 2011, their Fig. 2). However, up- 
dating the air-sea roughness parameterization led to 
only a 45% reduction of the original bias in surface wind 
speed. Specifically, the bias in zonally averaged South- 
ern Ocean surface winds in the CONTROL experiment 
is 2.75 m s -1 , but the reduction shown in Garfinkel et al. 
(2011) and discussed here is only 1.2 ms -1 . The sub- 
sequent stratospheric improvement is 25%-50% of the 
original stratospheric bias. 

In this section, we demonstrate that if the surface wind 
speed bias was removed entirely, 50%-100% of the 
stratospheric biases would be alleviated. We have per- 
formed a 27-yr sensitivity model experiment where the 
drag at the air-sea interface is increased far beyond what 
is observationally or theoretically justified, such that the 
surface wind speeds become realistic over the Southern 
Ocean. This is achieved by doubling the drag coefficient 
[C D in Eq. (2) of Garfinkel et al. (2011)] from 45°S and 
poleward after MOST has been applied, rather than 
updating the air-sea roughness coefficients as in Table 2 


of Garfinkel et al. (2011). Figure 13 shows the strato- 
spheric seasonal evolution in these experiments. The bias 
in column ozone in SON is removed, while the mid- 
winter vortex is unchanged. Overall, the decrease in the 
model biases in the stratosphere are approximately 
double that in the more realistic experiments shown in 



-60 


60 120 180 240 300 

longitude 


Fig. 11. Averaged 250-hPa height (m; Fig. lOd), vorticity (s -1 ) 
averaged from 150 to 200 hPa, and the total (azonal + zonal) 
vorticity forcing (s -2 ) averaged from 150 to 200 hPa (both the 
transient eddies themselves and the interaction of transient eddies 
with the mean flow) for UPDATED minus CONTROL. Each 
quantity has been averaged between 50° and 70°S. Heavy smooth- 
ing has been applied, though results are qualitatively similar with- 
out smoothing. 
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Fig. 12. Schematic of the response to increased surface rough- 
ness over the Southern Ocean. When the roughness is updated, 
drag is increased in a nearly zonally symmetric manner, as in- 
dicated by the orange easterly arrows in both the Indian Ocean and 
Pacific sectors. However, eddies (represented by dark gray arrows) 
act to mitigate this easterly anomaly. Because eddies are stronger 
climatologically, and also in response to an external forcing, in the 
Indian Ocean sector, the dark blue arrow is longer in the Indian 
Ocean sector than in the Pacific sector. The net effect is that the 
deceleration of the wind, which is represented with a purple-gray 
arrow, is larger in the Pacific sector than in the Indian Ocean sector. 

Fig. 3, so that there appears to be a linear relationship 
between improving the surface wind speed bias and the 
stratospheric SON biases. 

These experiments show that in a sensitivity (i.e., 
supertuned) experiment in which the surface winds are 
forced to become realistic, the improvement in the 


Polar cap temperature, lOOhPa 



stratospheric circulation is even more pronounced than 
when the available air-sea roughness data are used to 
update the model. These experiments suggest that ad- 
ditional, observationally justified, improvements to the 
model’s air-sea drag parameterization that result in 
further improvements to the surface wind speed may 
also result in further improvements to the strato- 
spheric flow. 

7. Discussion and conclusions 

An update to the air-sea roughness parameterization 
in NASA’s Goddard Earth Observing System Chemistry- 
Climate Model (GEOSCCM) leads to a decrease in 
model biases in both the troposphere and stratosphere. 
More than half of the model’s SH stratospheric biases 
appear to be related to air-sea roughness and surface 
wind speed biases. Increased surface roughness leads to 
an improved stationary wave pattern and enhanced up- 
ward wave activity flux entering the stratosphere in 
September and October. The increased wave driving 
leads to a more realistic Southern Hemisphere vortex: the 
vortex is no weaker in midwinter yet breaks up earlier in 


(b) vortex breakup (U<0) at 60S 



Polar cap ozone inter-annual variability 



Fig. 13. As in Fig. 3, but for the experiment in which the surface drag is increased beyond that which is justified by the 
available data on the air-sea interface, but in which the surface wind speed bias is entirely removed. 
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spring, which leads to improved polar cap temperatures 
and ozone from November through January. 

A possible mechanism whereby increased surface 
roughness leads to a zonally asymmetric stationary wave 
response is presented. Barnes and Garfinkel (2012) 
showed that increased drag on the mean flow leads to 
increased eddy momentum flux convergence to the jet 
core. Here, we demonstrate that this effect occurs pri- 
marily in regions where eddies are strong climatologi- 
cally (i.e., in the Indian Ocean sector and not the Pacific 
sector of the Southern Hemisphere). The wave-1 pattern 
of the eddy momentum flux convergence anomaly leads 
to a wave-1 pattern of the reduction in zonal wind and, 
by geostrophy, to a wave-1 geopotential height response. 
As the wave-1 anomaly is in phase with the climatological 
tropospheric stationary waves, wave-1 heat flux in the 
troposphere and the stratosphere is increased. Future 
work with simpler models is necessary to confirm this 
mechanism. 

Other studies have examined the impact of changes in 
other parameterizations on vortex biases. For example, 
McLandress et al. (2012) find that by parameterizing 
the orographic gravity waves associated with subgrid 
scale islands in the Southern Ocean, the delayed vortex 
breakup in the Canadian Middle Atmosphere Model 
(CMAM) is improved. However, they note that it is not 
possible to simultaneously hasten the springtime breakup 
without weakening the midwinter vortex. In contrast, 
we find that this change in air-sea roughness weakens 
the vortex in spring and early summer only. As resolved 
waves tend to only weakly influence the vortex in mid- 
winter (Taguchi and Yoden 2002, and references therein), 
it is perhaps to be expected that an improvement in- 
volving resolved waves (i.e., air-sea roughness) would 
be more effective in alleviating the late-breakup bias 
without exacerbating the midwinter bias. Furthermore, 
Hurwitz et al. (2010) attribute the delayed breakup in 
GEOSCCM to overly weak SH lower-stratospheric 
heat flux, and in particular, to weak stationary planetary 
waves in the troposphere, and here we have shown that 
the updated roughness parameterization leads to im- 
proved stationary waves in the midtroposphere (cf. 
Figs. 6 and 9). Finally, the newest implementation of 
the GEOSCCM includes orographic gravity wave drag 
associated with isolated small mountains in the subpolar 
SH (Molod et al. 2012, section 2.4), but the model 
biases persist. While an unresolved wave source might 
be important in other models (e.g., CMAM), an im- 
proved representation of resolved planetary waves in 
both the troposphere and stratosphere appears to be 
more important for GEOSCCM. 

The sensitivity of the stratospheric circulation to 
the air-sea roughness parameterization highlights the 


complexities in the atmospheric system. While it is dif- 
ficult to conclusively isolate our mechanism linking air- 
sea roughness to stationary waves in a comprehensive 
chemistry-climate model, it is remarkable that a rela- 
tively simple change in the surface-layer parameteriza- 
tion leads to such an improvement in the stratospheric 
circulation. As many other CCMs use similar schemes 
for their air-sea roughness parameterization and have 
similar biases in the stratosphere (and also in Southern 
Ocean winds; Son et al. 2010, their Fig. 10), the results 
from GEOSCCM may be relevant to other models as 
well. 
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